Electrical properties of improper ferroelectrics from first principles 
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We study the interplay of structural and polar distortions in hexagonal YMnC>3 and short-period 
PbTi03/SrTiC>3 superlattices by means of first-principles calculations at constrained electric dis- 
placement field D. We find that in YMnC>3 the tilts of the oxygen polyhedra produce a robustly 
polar ground state, which persists at any choice of the electrical boundary conditions. Conversely, 
in PTO/STO the antiferrodistortive instabilities alone do not break inversion symmetry, and open- 
circuit bundary conditions restore a non-polar state. We suggest that this qualitative difference 
naturally provides a route to rationalizing the concept of "improper ferroelectricity" from the point 
of view of first-principles theory. We discuss the implications of our arguments for the design of 
novel multiferroic materials with enhanced functionalities, and for the symmetry analysis of the 
phase transitions. 

PACS numbers: 71.15.-m,77.84.-s,77.55.Nv 



I. INTRODUCTION 

Multifunctional oxide materials, characterized by coex- 
isting ferroic orders of different physical origin, are cur- 
rently the object of intense exploration. Particularly en- 
ticing for their potential technological applications and 
fundamental interest are the so-called multiferroic mate- 
rials, where spontaneous magnetic and ferroelectric or- 
ders coexist in the same phase. Unfortunately, few mul- 
tiferroic oxides are known, and their magnetoelectric 
coupling is usually too small (or occurs at exceedingly 
lower temperatures) for practical purposes. This is typ- 
ically explained as follows: in the most popular ferro- 
electric oxides, such as BaTiOs, the spontaneous polar- 
ization emerges from the hybridization between occupied 
0(2p) orbitals and empty 3c? orbital of the transition- 
metal cation. This <i°-ness condition typically rules out 
magnetism. Although exceptions have been identified^, 
ferroelectricity and magnetism often appear to be con- 
traindicated. 3 

One way to achieve a strong coupling between electri- 
cal and magnetic degrees of freedom, is to resort to a 
more restrictive class of materials, called "improper fer- 
roelectrics" , where the polarization is induced by com- 
pletely different microscopic mechanisms. 4 (Another less 
restrictive mechanism involves the idea of phase com- 
petition^^.) Strictly speaking, the distinction between 
proper and improper ferroelectrics is made on symme- 
try grounds, and is based on the dynamics of the phase 
transition from a high-temperature centrosymmetric to a 
lower-temperature polar phase. In proper ferroelectrics 
(PF), such as PbTiC>3 or BaTi03, the polar distortion 
acts as the primary order parameter in the phase transi- 
tion; conversely, in improper ferroelectrics (IF) the polar- 
ization P is slave to another, non-polar, order parameter 
of different physical origin through an energy term linear 
in P. Depending on the material these can be the rota- 
tion/tilts of the oxygen polyhedra network, or the mag- 



netic order associated, e.g. with spin cycloids. Unlike 
covalency-driven polar distortions, these alternate mech- 
anisms can potentially lead to strong magnetoelectric 
couplings. For example spin-driven improper ferroelec- 
tricity in TbMnOa 7 naturally couples magnetic and polar 
degrees of freedom, resulting in a remarkable magneto- 
electric effect, while structurally-driven improper ferro- 
electricity in YMnO^ leads to an unusual locking of po- 
lar and antiphase structural domain walls^. These recent 
discoveries suggest an opportunity to achieve entirely 
new classes of strongly coupled multifunctional materials 
with novel properties based on improper ferroelectricity. 
As we said above, however, improper ferroelectrics are 
not abundant in nature. As a result material design rules 
have remained elusive and the discovery of new multifer- 
roics occurs serendipitously. 

This situation is now rapidly changing thanks to the 
increasing popularity of density-functional theory calcu- 
lations, which have led the quest for new oxide materi- 
als in the past few years. An important breakthrough 
in the fielcU^ occurred by the recent discovery of a new 
type of improper ferroelectricity in artificially layered 
PbTi0 3 /SrTi0 3 (PTO/STO) superlattices. This sys- 
tem displays an unusual form of lattice-driven instabil- 
ity, involving the interplay of P with octahedra rotations 
and tilts. More specifically, the free energy of this sys- 
tem was shown to display an unusual trilinear coupling, 
^$i$2-P, between the polarization and two non-polar 
antiferrodistortive (AFD) modes of different symmetry, 
$i and <I>2. This work has paved the way towards the 
rational design of lattice-driven improper ferroelectric^ 
with a large spontaneous P and has lead to the discov- 
ery of a novel class of (CaMnOs^CaO strongly coupled 
multiferroics 1 - 2 -. Given the versatility of perovskites, this 
novel class of IFs are ripe systems to explore for new 
multifunctional phenomena. 

There are however important questions that still need 
to be answered from the point of view of fundamental 
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theory. For example, the origin of the polarization is 
not always clear if in addition to the two AFD instabil- 
ities - which from the trilinear coupling could lead to 
a P ~ $i$2 _ a proper ferroelectric instability exists. 
This issue is ultimately tied to the question of which two 
order parameters can be identified as inducing the dis- 
torted state (e.g., in this case $1 ~ $2-P is an equally 
valid statement). This question of primary order pa- 
rameters determines to what extent the electrical and 
dielectric properties behave like conventional IFs. This 
is of central importance for the applicability of multi- 
functional materials in actual devices. Unfortunately the 
ability to clearly address this using symmetry arguments 
and/or conventional first-principles calculations alone is 
difficult. What is clearly missing is a generally applicable 
and reliable "test" within a zero-temperature calculation, 
which would unambiguously identify whether the prop- 
erties of this class of (hybrid) IFs behaved more like a 
conventional IF or more like a proper ferroelectric. 

Here we compare from first principles the electrical 
properties of hexagonal YMn0 3 to that of the PTO/STO 
superlattice. We illustrate the key fingerprints of im- 
proper ferroelectricity by using a recent advance in first- 
principles methods^ that treats the electrical displace- 
ment, D, as the fundamental variable. This allows us to 
directly compute the spontaneous polarization in open- 
circuit, PoCi within a periodic bulk calculation, avoiding 
the need for cumbersome slab geometries (where the re- 
sults inevitably depend on the details of the surfaces; a 
common source of artifacts in the calculation). We iden- 
tify the presence of a nonzero and switchable Poc a s 
a necessary (although generally not sufficient) condition 
that must be satisfied by a material in order for it to 
be classified as improper ferroelectric. This property is 
enjoyed by YMnC>3 (albeit with a substantially smaller 
spontaneous Poc than reported previously), but not by 
the 1/1 PbTiC>3/SrTi03 super lattices considered here. 
We rationalize this unexpected result in terms of higher 
order couplings between the AFD mode(s) and the polar- 
ization. We briefly discuss the implications of these result 
for the future design of novel improper ferroelectric ma- 
terials. This approach puts the concept of improper fer- 
roelectricity on firmer conceptual grounds from the point 
of view of ab-initio theory, providing a useful tool for the 
rational design of new "hybrid"— IFs. 



II. COMPUTATIONAL METHODS 

Our calculations are performed within the local- 
density approximation (LDA) of density-functional the- 
ory and the projector-augmented- wave (PAW) method 14 , 
with a plane- wave cutoff of 60 Ry. For YMnC>3 , we ap- 
ply a Hubbard correction (LDA+U) of U=6 eV to the 
Mn 3c? states, in order to obtain an insulating ground 
state. We further assume a collinear antiferromagnetic 
configuration of the spins, and use a 4x4x2 Monkhorst- 
Pack 15 (MP) sampling of the 30-atom hexagonal unit 
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FIG. 1: Electrical properties of YM11O3. The energy (a), 
internal electric field (b), K3 mode amplitude (c) and out-of- 
plane strain (d) are plotted as a function of the macroscopic 
D field. Symbols refer to first-principles data; the continuous 
lines are SPLINE fits to the symbols. 



cell to perform the Brillouin-zone integrals. We use stan- 
dard non-spin polarized LDA for PTO/STO, together 
with a 4x4x3 MP grid to sample the Brillouin zone of 
the x x 2 direct-space cell (20 atoms total, dou- 
bled in-plane to accommodate the antiferrodistortive ro- 
tations of the oxygen octahedra). Cell relaxations are 
performed by using the calculated stress tensor, which 
we correct for the Pulay error by applying a constant 
negative pressure. For YMn03, first we optimize the 
structure of the ferroelectric phase, and then fix the in- 
plane lattice parameter to the equilibrium value in zero 
applied field. For PTO/STO we fix the in-plane lattice 
parameter to the equilibrium theoretical value of cubic 
SrTiOs, as it has became common practice in simulations 
of layered perovskites. (We relax only the out-of-plane 
lattice parameter in response to the electrical perturba- 
tions, as typically these materials are grown as epitaxial 
films on a thick substrate, which prevents their in-plane 
relaxation.) The macroscopic electrical degrees of free- 
dom are dealt with by using the methods described in 
Ref. [lU In all cases we enforce rather stringent thresh- 
old for structure relaxations, as the polarization and the 
internal electric fields are both very sensitive to minimal 
displacements of the ions. More specifically, we calculate 
the ground state of both materials as a function of the 
macroscopic electric displacement field, and analyze the 
results in terms of energy, internal electric field, ampli- 
tude of the relevant structural (zone-boundary) modes 
and out-of-plane (electrically induced) strain. 
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III. RESULTS 
A. YMn0 3 

In a proper ferroelectric, open-circuit boundary condi- 
tions suppresses single-domain ferroelectricity and drives 
domain formation. Imposing open-circuit boundary con- 
ditions corresponds to adding to the free energy a positive 
term that is quadratic in P— , 

&E=^(D-AkP)\ (1) 

where f2 is the volume and D is the electric displacement 
of the high-temperature paraelectric phase, which can be 
taken as zero. Remarkably, Sai et alJ£ demonstrated by 
means of first-principles calculations that a monodomain 
state with non-zero P can persist in YMnC>3 even in the 
extreme limit of an unsupported film without electrodes. 
Because the polarization couples linearly to the primary 
order parameter in an improper ferroelectric^, any IF 
must enjoy such a property. 

So are there clearly identifiable features in the elec- 
trical diagram that distinguish an improper ferroelectric 
from that of a proper ferroelectric? In FigureQJa) we plot 
the energy as a function of the electric displacement field 
for YMnC>3. As expected from a ferroelectric compound, 
there are two symmetric potential wells, whose minima 
correspond to a spontaneous polarization of P s = 7.1 
/LtC/cm 2 . This is slightly larger than previously reported 
values 8 -*^; we ascribe the small difference to the different 
pseudopotentials used, and to the different (real-space) 
approach used here to compute the macroscopic polar- 
ization 17 . Remarkably, the electrical equation of state 
is a multi-valued function in the region near the origin. 
Moreover, the curves belonging to the two symmetrical 
wells cross each other in D — 0, rather than converging 
smoothly into a saddle point. Both features are in stark 
contrast to the known behavior of "standard" ferroelec- 
tric materials that display a proper ferroelectric transi- 
tion such as PbTi0 3 or BaTi0 3 . 

This is not completely unexpected, as the physical 
origin of ferroelectricity in YMn03 differs from that of 
BaTiC>3. YMn03 is a classic improper ferroelectric, in 
which the spontaneous P is driven by an unstable, non- 
polar "trimerization" (A3) mode such that P oc Af. 
Thus, rather than destabilizing the polarization, the A3 
distortion pushes the zone-center polar modes away from 
their centrosymmetric equilibrium configuration, gener- 
ating a polarization in the process. In other words, the 
A'3 mode acts like a switchable "geometric field"— that, 
depending on its orientation, produces two inequivalent 
(but related by a mirror symmetry operation) states of 
the system with opposite tendency towards a polar dis- 
tortion. Either state shows an almost perfectly linear di- 
electric behavior [see FigureQJb), the electric field follows 
an almost perfect straight line], with an average dielec- 
tric constant €33 ~ 13, and an extremely weak piezoelec- 
tric response, of the order of CZ33 ~ 1 pm/V. Note that 
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TABLE I: Properties of the relevant phases in PTO/STO. 
The spontaneous polarization P z is in /xC/cm 2 ; the angles <j> 
are in degrees; energies AE (referred to the PA/mmm phase) 
are in meV. 



the amplitude of the K3 phonon mode does not change 
much across the whole range of D values, max ±5 % with 
respect to the ferroelectric ground state. This is a conse- 
quence of the relatively weak coupling (only third-order) 
between the K3 mode and P, and to the comparatively 
high energy associated to the K 3 structural transition, 
which is more than three times higher than the highest 
electrostatic energy, about 0.15 eV, obtained in our sim- 
ulations. In passing, we note that a similar kink in the 
electrical equation of state at P = was previously re- 
ported for another ferroelectric material, KNO3, where 
there is also a structural mode that produces a bistable 
state in the origi n 18 ' 19 . We identify this characteristic 
kink, which is a consequence of the above-mentioned "ge- 
ometric field" , as the ab-initio fingerprint of improper 
ferroelectricity. 

Our calculations indeed show, in agreement with 
Ref. [l6| that there is a non-zero polarization in open 
circuit, Poo, which corresponds to the D = point 
in the electrical diagrams of Figure [1] Since D = 0, 
Poc — —Eoc/^i where £oc — —0.62 GV/m is the 
internal field; this yields Poc = ±0.55 ^tC/cm 2 . (The 
ratio between the spontaneous polarizations in short cir- 
cuit and open circuit corresponds to the average dielec- 
tric constant of YMn03 within this region of the electri- 
cal diagram: P s /Poc — £33 ~13.) Note that these val- 
ues depart significantly from those (Pop =6.1 /iC/cm 2 , 
£oc = 0.25 GV/m) reported in Ref. fig. To account for 
this discrepancy, we can speculate that the unsupported 
films of Ref. [l6| are not under ideal open-circuit bound- 
ary conditions: The metallic oxygen-terminated surfaces 
would act effectively as electrodes, that partially screen 
the polarization charges of YMn03. This highlights the 
importance of the analysis introduced in this work for 
understanding the intrinsic properties of improper fer- 
roelectrics. We shall come back to this point in Sec- 
tion [TTTCl where we show how one can predict, from the 
calculated bulk data alone, the behavior of a film under 
any type of electrical boundary conditions, ranging from 
open-circuit to close-circuit, and including all intermedi- 
ate "imperfect screening"— regimes. 
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FIG. 2: Electrical properties of the PTO/STO superlattice. 
The energy (a), internal electric field (b), AFD mode ampli- 
tudes (c) and out-of-plane strain (d) are plotted as a function 
of the macroscopic D field. 



B. PbTi0 3 /SrTi0 3 superlattices 

We now consider the PTO/STO superlattice. In 
the high-symmetry P4/mmm reference phase we find 
three unstable distortion modes, in agreement with 
Ref. [Tot a zone-center polar mode, plus two symmetry- 
inequivalent AFD rotations of the oxygen octahedrals 
along the z axis, either out of phase (AFD zo ) or in 
phase (AFD z i). (Note that we restrict our analysis to 
modes with rotation/displacement axis along z.) By 
freezing in each of these instabilities separately and then 
performing full structural relaxations we obtain three 
new phases of lower symmetry labeled FE 2 , AFD 2 j, 
and AFD 2Q in Tabic |U Furthermore, we obtain a fourth 
phase, (AFD/FE) 2 , in which all three instabilities coex- 
ist. While we find some discrepancies at the quantitative 
level between our results and that of Ref. the level 
of qualitative agreement confirms their main conclusion. 
The energy gains associated to ferroelectric transitions 
[PA/mmm ->• FE 2 , AFD zo — >(AFD/FE) 2 ] are neverthe- 
less systematically larger by a factor of 3-4 here (respec- 
tively -3 and -8 meV in Ref. [13, to be compared with 
-12 and -56+32=-24 meV in Table HJ}, while the AFD zi 
and AFD 2Q instabilities appear systematically less pro- 
nounced (-28 and -73 meV in Ref. Q3). Ferroelectric 
perovskites are known to be especially sensitive to the 
computational approximations and these differences can 
be traced back in the different pscudopotentials used in 
both calculations and in the subsequent different in-plane 
lattice constant at which the calculations have been per- 
formed. 

The large energy gains associated with the ferroelectric 
transition that we calculate suggests in fact that the po- 
larization might be one of the primary order parameters 
(the other being AFD zo ) rather than the two rotation 




FIG. 3: Two-dimensional energy landscape of the AFD modes 
in the PTO/STO superlattice around the P4/mmm reference 
structure. Horizontal {4> zo ) and vertical axes are the 

mode amplitudes in degrees. The unit of the color bar is 
meV. 



distortions. Applying the approach outlined in this work 
should clarify this point. In Figure [2] we summarize the 
non-linear electrical and structural response of the super- 
lattice as a function of D, analogous plots as in Figured] 
Unlike YMn0 3 , PTO/STO appears to behave in all re- 
spects like a standard ferroelectric material. The double- 
well potential in Figure [5Ja) is a single- valued function, 
diffcrcntiable at the saddle point, qualitatively similar 
to typical soft-mode driven systems like PbTi03 and 
BaTiOa. This was to some extent unexpected : recent 
reports of improper ferroelectricity on this same system 
suggested that its properties might closely mirror those of 
YMnOa 10 . To gain insight into our apparently contradic- 
tory results, it is useful to look at the response of the AFD 
modes - which interact nontrivially with the ferroelectric 
polarization - to an applied Z?-ficld. In Figure [2][c) we 
show the amplitude of the out-of-phase AFD zo and of 
the in-phase AFD z i modes as a function of D. At D = 0, 
only a pure AFD zo distortion is present (i.e., this config- 
uration is identical to the pure AFD 2Q phase reported in 
Table H|, while the AFD Z ; distortion and the polarization 
are absent (AFD 2G is a centrosymmetric space group.) 
This explains why we do not observe a linear behavior 
of the energy around D = in Figure (Ha) . Only when 
we move away from D — do we see a linear increase of 
AFD Z ; , the signature of the tri-linear coupling discussed 
in Ref. O 

To illustrate why the AFD zo and AFD Z ; do not co- 
exist at D = 0, we performed a series of calculation 
where we freeze in various amplitudes of either mode, 
starting from the relaxed superlattice in the reference 
Pi I mmm symmetry. In Figure [3] we plot the energy as 
a function of the two independent amplitudes (4> zo , <p z i) 
within the restricted two-dimensional subspace spanned 
by these modes. As speculated above, both modes are 
unstable in the reference structure, but the condensation 
of either mode quickly suppresses the instability with re- 
spect to the other (consistent with the results of RefJ^) . 
The global minima of the two-dimensional energy land- 
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Inverse interface capacitance density (m /F) 

FIG. 4: Equilibrium D for a thin-film PTO/STO (solid line) 
or YMn03 (dashed line) capacitor as a function of X e s/N or 
Ci 1 /N. The shaded area refers to a hypothetical overscreen- 
ing regime with A cr r < 0. The dotted line is a guide to the 
eye, to highlight the linear behavior of D in PTO/STO for 
X e f[/N < 0.025 A. The inset shows the absence of "critical 
thickness" in YMn03 in the limit of large X e s/N. 

scape corresponds to pure AFD zo states (the out-of-phase 
instability is much stronger than the in-phase one), in 
agreement with the D = ground state found in our 
electrical diagram of Figure^ a). 

In the following Section we further substantiate these 
points by comparing the performance of hypothetical 
PTO/STO- or YMn0 3 -based capacitors with real elec- 
trodes (i.e. under imperfect screening conditions). 



dU jdD = (recall that the potential V is related to U 
by SV(D) = dU/dlM), 

on 

V(D) = NV B uik(D) + — = 0. (3) 

This can be written more compactly in terms of the in- 
terfacial "effective screening length" A c ff = eo/Ci as 

Acft _ ^Bulk(-P) /js 

e Q N ~~ 2D U 

By inverting Eq. U we can predict the behavior of a thin 
film of arbitrary thickness ./V and with an arbitrary A e g- 
(the ground state depends on X e s/N only). 

In Fig. SJa) we plot the resulting equilibrium D of a 
hypothetical N = 1 capacitor as a function of \ e g. (We 
also consider small negative values of X e s/N to account 
for a possible overscreenin g 21 ' 23 .) PTO/STO (solid line) 
shows the typical destabilization of the ferroelectric state 
(critical thickness for ferroelectricity) under an imperfect 
screening regime; this occurs at A c ff /iV=0.047 A. Con- 
versely, YMn03 appears almost unsensitive to the electri- 
cal boundary conditions, with a spontaneous polarization 
that persists up to arbitrarily large values of X e s/N. The 
limit X c g/N —¥ oo corresponds to open-circuit bound- 
ary conditions, where D tends to zero asymptotically as 
~AT/A eff . 



C. From bulk to thin- film capacitors 

So far we have only discussed the bulk electrical prop- 
erties of YMn0 3 and short-period PTO/STO superlat- 
tices. However, for practical applications in thin-film ca- 
pacitors, it is necessary to take possible finite size effects 
into account. In particular, the imperfect screening at 
a realistic film/electrode interface is known to affect the 
stability of the ferroelectric state by altering the electri- 
cal boundary conditions of the fil m 20 ' 21 . In the following 
we shall demonstrate that the impact of such interface 
effects can be directly predicted from the bulk equations 
of state of Fig. [1] and [2j i.e. without doing explicitly a 
demanding calculation of an actual capacitor. 

We assume two symmetric electrode/ferroelectric in- 
terfaces, which behave as two thin layers of linear di- 
electric interposed between the insulator and the metal. 
These can be characterized indifferently by an interface 
capacitance density per unit area Ct, or by an effective 
screening length A e fi = e /Cj^. Then, the electrical 
equation of state of such a device can be written as^ 

U(D) = NU Bn]k (D)+2—, (2) 

where Usuik is the energy as in Fig. QJa) and [U[a); S 
is the cell cross-section and N is the thickness of the 
film expressed in number of bulk units. The mini- 
mum of U(D) within short-circuit boundary conditions 
(zero applied potential) is given by the stationary point 



D. From first principles to Landau theory 

While our first-principles results yield a fairly detailed 
picture on the electrical behavior of the compounds un- 
der consideration, it is very useful to recast these data 
in terms of a simplified physical model, where the role 
played by the active degrees of freedom emerges clearly. 
Also, most researchers in the ab-initio community are 
used to analyzing the behavior of ferroelectrics by per- 
forming calculations under zero external field, where the 
energy landscape as function of the active structural dis- 
tortions is explored by freezing in various amplitudes of 
selected imaginary-frequency phonons. In this Section 
we shall trace a formal link between the constrained- D 
method used in this work and the more familiar under- 
standing of ferroelectrics, within either Landau theory or 
the frozen-phonon approach. 

Consider the zero-field first-principles Kohn-Sham to- 
tal energy, written as a function of the electronic degrees 
of freedom ip and the ionic positions u (we omit the ex- 
plicit discussion of the strain to avoid complicating the 
notation), -Eks({"0}j i u })- Around the high-symmetry 
configuration, one can construct a harmonic expansion 
in terms of the N lattice eigenmode amplitudes £j , 

£W{0)~ai£i 2 + --- + a^, ( 5 ) 

with the assumption that all vanish in the reference 
structure. The polarity of the phonons can be expressed 
via the mode effective charge, Zi = fl dP/d^i (we assume 
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transversal zone-center modes, and we restrict the anal- 
ysis to one component of the polarization only) . Assume 
that the lattice distortions do not significantly affect the 
purely electronic dielectric constant, e^, which we con- 
sider linear in the applied external field; then we can 
readily write the constrained- 1? functional 

[/({£}, D) = E KS ({0) + ^— — (D - P 5 ) 2 , (6) 

where is the lattice-mediated contribution to the po- 
larization, 

P 6 = ^iZi + ...£nZ n ). (7) 

The ground state at a given D is then obtained by tak- 
ing the global minimum with respect to all degrees of 
freedom, 

U(D) =min U({Q,D) (8) 

Now, to make progress from here it is useful to note 
that not all phonon modes are equally important to de- 
scribe the material. Most modes are hard, with little 
polar activity. Only few low-energy phonons display a 
marked anharmonic behavior or are unstable; these are 
the modes that we wish to focus our attention on. Typ- 
ically, the hard modes are discarded from the Hamilto- 
nian, or better, are incorporated into an effective back- 
ground constant; conversely, the low-energy modes are 
kept explicit, and higher-order terms are introduced in 
the Hamiltonian to ensure their proper description. As- 
sume for example that only one mode, £ s , is "soft", and 
that it can be reasonably well described by a quadratic 
plus a quartic term, 

£Ks(6) = a s £ s 2 + /^ s 4 . (9) 

Assume also that the phonon-phonon interaction is small, 
and that the electronic susceptibility is not significantly 
affected by £ s . Then the constrained- D functional can 
be simplified by taking the minimum with respect to the 
hard modes only, 

Ufo,D) = mm Um,D) 

= a s e s +ne s +^-{D-p s f. (to) 

Here we have introduced the background dielectric con- 
stant^ cb, which is given in terms of e^, the (hard) 
mode charges Zi and their eigenfrequencies u>i, 

Open-circuit boundary conditions correspond to setting 
D = in the above equation. In such a regime it is easy 
to see that the electrostatic energy term in U is just an 



additional quadratic (and positive) term in the soft mode 
amplitude, of the type A£ 2 . Typically, this additional 
term completely suppresses the ferroelectric instability: 
a s + A > 0. 

The above discussion establishes an important result: 
if we know the approximate low-energy Hamiltonian in 
zero field, the "background permittivity" and the effec- 
tive charges of the active modes (all these ingredients can 
be indifferently extracted from ab initio calculations or 
from a phenomenological analysis) we can immediately 
calculate the ground state at a specific value of D. In 
other words, the so-called constrained- 1? method is by 
no means specific to ab initio calculations, and can be 
applied just as well to effective Hamiltonians and phe- 
nomenological descriptions. Indeed, mapping the proper- 
ties of a given material in "electric displacement space" 25 
may be a powerful tool to validate second-principles or 
macroscopic approaches against the microscopic first- 
principles theory - if the descriptions are of the same 
quality, the results should coincide. In the following 
we discuss how these ideas apply to our two testcases, 
YMn0 3 and PTO/STO. 

In YMn03 the low-energy Hamiltonian can be written 
as a function of a polar (but not soft) T mode and by an 
unstable if-point mode^ 

+ Pkk& + Ptk&& + 7fr& (12) 

The electrical variables have little impact on the struc- 
tural if -point mode, that acts like a "switch" (see Fig.[TJ), 
and can therefore be approximated by a constant with ± 
sign, 

£«s(fr)~ar£r+&r±7fr. (13) 

Typical values of the parameters have been reported by 
Fennie and Rabe — . The biquadratic coupling renormal- 
izes the T-mode frequency by further hardening it - the 
new quadratic coefficient is given by a — ar + /3. The 
last term is responsible for the geometric field. If we plug 
this Hamiltonian into the electrostatic energy functional 
we obtain 

U(^,D) = a^±^ + -^-(D-Pr) 2 , (14) 

2eo£B 

where Pr = Zp£,r/Q- This is a second-order polynomial 
that can be solved analytically for the ground state at a 
given D, 

U(D) = -^—(D±P ) 2 (15) 

^C0 e TOT 

where €tot is the total dielectric constant, including the 
contribution from the T mode, and P$ is the sponta- 
neous polarization of YMJ1O3 in zero field. This equation 
accurately describes the results shown in Fig. [I] where 
YMn03 emerges from the calculations as a linear di- 
electric with £tot ~ 13 and a switchable polarization 
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Po = 7.1 C/m 2 . In particular, if we impose open-circuit 
boundary conditions by setting D = 0, the electric field 
is 



dU 

£oc = dD 



D=0 CoETOT 



(16) 



and the polarization is Poc = ~ e o^oc — T-Po/ e TOT) 
consistent with our analysis of the first-principles data. 

The case of PTO/STO is slightly more complicated, 
due to the presence of three active modes instead of two, 
but we can follow a similar line of reasoning. The ba- 
sic physics can be understood in terms of three phonon 
modes of the high-symmetry PA/mmm reference struc- 
ture: a polar T mode £r, and two zone-boundary struc- 
tural modes that correspond to out-of-phase and in-phase 
rotations of the oxygen octahedra with respect to the z 
axis, (j) zo and <p Z i. The energy in zero field can be ex- 
panded as 



i) 



A 

A . 

BqT 



2 

oVzo 
2 



+ B CI 

+ Bi,,,_, 



5 rr?r + 

'oo'Pzo 
'u^zi + 



B n 



(17) 



The first three lines correspond to the double- well poten- 
tial associated with each individual mode (or, a zo , a Z i < 
0); the fourth line provides the lowest-order biquadratic 
couplings responsible for the pair-wise competition or co- 
operation between the modes; the last line is the tri-linear 
coupling, that was previously identified as key to a novel 
type of improper ferroelectricity in this system. As the 
rotations do not directly contribute to P (their dynamical 
charge is zero), the simplified constrained- D functional 
can be written as 



U(£r,4>zo,(t>zi,D) = -^KS (£r, </>2o, 0zi) + 



2eo£B 



{D-Pvf. 
(18) 

At D — 0, the electrostatic energy term reduces to 7£p 
with 7 a positive constant; this effectively renormalizes 
the Ay coefficient to Ay = Ar+7. Recall that Ay is nega- 
tive - in our ab initio model PTO/STO has a ferroelectric 
instability even in absence of rotations. Conversely, Ay 
is positive and large. Indeed, the electrostatic term in U 
induce a shift in the frequency of the zone-center optical 
modes which become longitudinal rather than transverse, 
and the LO-TO splitting in perovskite titanates is enor- 
mous^ As Ay is positive, the only term in the Hamil- 
tonian that can induce a polarization in open-circuit is 
precisely the trilinear one C(j)zo4>ziS,v- Furthermore, as 
Ay is large, £r will be virtually suppressed from the low- 
energy Hamiltonian, and only a simultaneous presence of 
zi and zo rotations (at £r = 0) will produce a geometric 
field. In order to verify this possibility, we need to study 
the energy as a function of the AFD modes amplitude, 
i.e. the energy landscape that is sketched in Fig. 3. 



A a 




Bu 


Bi 


-2.1097 


-0.6912 0.0368 


0.0363 


0.1372 



TABLE II: Parameters of the model described in the text. 
Values of the parameters are in eV. The independent variables 
<f> are in degrees. 



In absence of £r, the low-energy Hamiltonian reduces 



to 



Eks(4>zo, <t>zi) = A 4> 2 Z0 + Ai4> 2 zi + 



B n 



Bj, 



B 



loVzo^zi " 



(19) 



The values of the five parameters are reported in Ta- 
ble [TTJ the maximum deviations with the actual first- 
principles data are lower then 1.5 meV across the whole 
two-dimensional energy landscape. 

The necessary and sufficient condition for having four 
minima where <p zo and <p Z i are nonzero is that the "pure" 
minima (where only one degree of freedom is active at 
the time) are in fact saddle points of the two-dimensional 
function E(<j> zo ,<f> zi ), 



A - 
Ai - 



AiBio 
2Bu 

A B{ 
^B nn 



< 

< 



(20) 
(21) 



It is easy to see that the first of these conditions is sat- 
isfied, i.e. the AFD Z ; configuration is not a true min- 
imum but a saddle point in this reduced configuration 
space. Conversely, the second condition is violated, i.e. 
the pure AFD zo state is indeed a global minimum, where 
the original in-phase AFD instability is suppressed. 



IV. DISCUSSION 

In the previous section we have clarified the physical 
ingredients that are responsible for the calculated electri- 
cal behavior of YMn0 3 and PTO/STO. However, a sub- 
stantial limitation of our work so far has been the com- 
plete neglect of temperature effects. As the concept of 
"improper ferroelectric" is intimately bound to the ther- 
modynamics of phase transitions, it is useful to discuss 
here whether we can infer a number of likely scenarios 
based on our zero-temperature results. (Of course, our 
arguments will remain limited to a qualitative level, as 
presently we don't have access to finite-temperature sim- 
ulations.) This point is particularly important in order 
to reconcile our findings (that a true "geometric field"— is 
not pres ent in PTO/STO) with the experimental data of 
Ref. Ilfj . where the typical signatures of improper ferro- 
electricity emerge clearly. Given the relevance of this sys- 
tem in the development of novel multifunctional materi- 
als, we shall focus the following discussion on PTO/STO. 
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Experimentally, PTO/STO undergoes a phase transi- 
tion at Tc = 500 K from a high-temperature paraelectric 
state to a polar state where the spontaneous P that in- 
creases linearly with decreasing T. To rationalize this 
phase transition, several scenarios are possible, depend- 
ing on the respective behavior of the three active degrees 
of freedom. First, we can imagine a transition where 
the AFD modes act as primary order parameters in the 
transition, and the polarization emerges as a secondary 
effect; this is the classic picture of improper ferroelec- 
tricity. This, however, does not appear consistent with 
our simulations. In absence of polarization, the two AFD 
modes are mutually exclusive, and they don't break in- 
version symmetry alone. In other words, if we suppress 
P there is no "geometric field" in the system. Second, we 
can imagine a transition where £r is the primary order 
parameter, and the AFD distortions are secondary. This 
is, however, inconsistent with the experiments of Ref. ITol : 
in a second-order phase transition involving only P, one 
would expect P oc \JTc — T, which differs from the re- 
ported linear behavior. 

Interestingly, a third possibility was proposed in a re- 
cent work;2i where analogous Hamiltonians involving a 
trilinear coupling were studied in a phenomenological 
framework. The scenario described in Ref. is that of 
an "avalanche" first-order transition, where all three de- 
grees of freedom condense simultaneously without any of 
them becoming unstable. This is a situation where the 
usual assumptions of Landau theory, and in particular 
the distinction between primary and secondary order pa- 
rameters, break down: neither the polarization nor the 
AFD modes would act as either master or slave. This 
possibility can be simply understood by constructing a 
"reaction coordinate" Q and by imposing that all three 
degrees of freedom must be proportional to Q. Then the 
Hamiltonian reduces to the following expression 

E KS (Q) = aQ 2 +(3Q A + 1 Q 3 . (22) 

All the degrees of freedom are individually stable above 
the transition, so a is positive, and its magnitude will 
decrease with T. j3 is necessarily positive at all tempera- 
tures, or thermodynamic stability is lost. One can easily 
see that if 7 is different from zero a secondary minimum 
will appear for sufficiently low T (small a), and will drop 
below the main minimum at Q = before the Q = min- 
imum becomes a saddle point. Thus, in presence of this 
trilinear term in the Hamiltonian, the transition is very 
likely to be first-order and involve all degrees of freedom 
simultaneously. 

Strictly speaking, even this hypothesis does not con- 
stitute an ideal match to the experimental results. A 
first-order phase transition necessarily involves a jump 
in the polarization from zero to a finite value at Tc, 
while a continuous curve was reported experimentally. 
One can imagine, however, that the transition might be 
only weakly first-order, and the discontinuity might be 
too small to be detectable in the experimental setup of 
Ref. [l^- Another possibility is that the depolarizing field 



produced by the electrodes might have an impact on the 
measured P, or more generally that other degrees of free- 
dom not considered here (e.g. nanodomains) might play 
a role in the transition. A last alternative is that stan- 
dard DFT calculations do not properly describe the rela- 
tive strength of the different instabilities although these 
methods have been successfully applied to many cases so 
far. 

Resolving these issues goes beyond the scopes of the 
present work; what we want to stress here is that the 
electrical analysis used in this work constitutes a pow- 
erful tool to study complex ferroelectric systems where 
many degrees of freedom interact nontrivially with the 
polarization. 



V. CONCLUSIONS 

In conclusion, we have shown that constrained- D 
calculations provide a powerful method to understand 
the electrical behavior of multifunctional oxide com- 
pounds. As practical examples, we discussed two differ- 
ent material systems: A prototypical improper ferroelec- 
tric, YMn0 3 , and short-period PTO/STO superlattices. 
Concerning YMn03, our results have evidenced a charac- 
teristic "kink" in the electrical equation of state (energy 
U as a function of the electric displacement D), which 
is responsible for the remarkable insensitivity of this ma- 
terial to the electrical boundary conditions^ Concern- 
ing PTO/STO, our results indicate that the electrical 
properties depends crucially on the interaction between 
the two non-polar AFD modes: Only if the AFD modes 
are the dominant instabilities and they do not compete 
with each other we expect the electrical properties to 
be qualitatively similar to those of YMn03. While our 
PTO/STO model does not appear to satisfy either con- 
dition, recent results^ give strong indications for such a 
scenario to occur in a related material, Ca3Mn2 07. To 
illustrate the consequences of these findings in a realistic 
capacitor configuration, we showed how the U{D) curve 
can be readily converted into a P(\ e ff/N) curve, where 
X e f f is the effective screening length of the electrode in- 
terface and N is the thickness. 

As an outlook, Ref. [13 has introduced us to the re- 
markable possibility that two AFD rotation distortions 
can induce ferroelectricity via a trilinear coupling with 
the polarization. This idea holds much promise to real- 
ize novel multifunctional phenomena, and has led to the 
discovery of an entirely new class of "hybrid" improper 
ferroelectric^. While rationalizing these breakthroughs 
by putting them on firm theoretical grounds, our work 
also illustrates some interesting new aspects of these sys- 
tems that might be useful for future applications. For ex- 
ample, it appears clear from our data that, in contrast to 
YMn03, in PTO/STO the AFD modes are very strongly 
coupled to the polarization. This implies that the con- 
trol of octahedral rotations by means of an applied field is 
indeed possible, and the constrained- D method appears 
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as the ideal theoretical tool to investigate (and engineer) 
this functionality in a variety of oxide systems. 
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